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Abstract. In this paper we present a Boltzmann type price formation model, which is 
motivated by a parabolic free boundary model for the evolution of the prize presented by 
Lasry and Lions in 2007. We discuss the mathematical analysis of the Boltzmann type model 
and show that its solutions converge to solutions of the model by Lasry and Lions as the 
transaction rate tends to infinity. Furthermore we analyse the behaviour of the initial layer 
on the fast time scale and illustrate the price dynamics with various numerical experiments. 
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1. A Boltzmann type model for price formation 
Market microstructure analysis studies the trading mechanisms of assets in (financial) mar- 



kets. According to O'Hara ( 1998 ) markets have two principal functions - they provide liquidity 
and facilitate the price. The evolution of the price is influenced by the trading system and the 
nature of the players, and is hence an emergent phenomenon from microscopic interaction, 
that calls for advanced mathematical modelling and analysis. 

While there has been a lot of research in economics and statistics related to either the micro- or 
macrostructure of markets, few investigations dealt with systematically analysing the micro- 



macro transition. Recent work by Lasry and Lions (2007) revived the mathematical interest 



by introducing a price formation model that describes the evolution of price by a system of 
parabolic equations for the trader densities (as functions of the bid-ask price), with the agreed 
price entering as a free boundary. The authors motivated the model using mean field game 
theory, but the detailed microscopic origin remained unclear. 

In this paper we provide a simple agent based trade model with standard stochastic price fluc- 
tuations together with discrete trading events. By modelling trading events between vendors 
and buyers as kinetic collisions we obtain a Boltzmann-type model for the densities. Then we 
prove rigorously that in the limit of large trading frequencies, the proposed Boltzmann model 
converges to the Lasry and Lions free boundary problem. We also analyse other asymptotics 
beyond the scales that the free boundary model can describe. Hence we provide a basis for 
deriving macroscopic limits from the microscopic structure of trading events, which allows for 
various generalisations to make the model more realistic. 
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ON A BOLTZMANN TYPE PRICE FORMATION MODEL 



To set up the model we consider the price formation process of a certain good which is traded 
between two groups, namely a group of buyers and a group of vendors. The groups are 
described by the two positive densities / = f(x,t) and g = g(x,t), where x denotes the bid 
resp. ask price and t the time. If a buyer and a seller agree on a price, a transaction takes 
place and the buyer immediately becomes a seller and vice versa. The transaction fee is 
denoted by a positive constant a G M. As a consequence the actual price for a buyer is x + a, 
therefore he/she will try to resell the good for at least that price. The profit for the vendor 
is x — a, therefore he/she is not willing to pay a higher price in the next trading event. 
The situation above can be described by the following Boltzmann type price formation model: 



ft(x,t) 



a 



fxx(x, t) - kf(x, t)g(x, t) + kf(x + a, t)g(x + a, t) 
> 

-g X x(x, t) - kf(x, t)g(x, t) + kf(x - a, t)g(x - a, t). 



(1.1a) 
(1.1b) 



9t(x,t) = 
with initial data 

f(x, 0) = //(ar) > 0, g(x, 0) = 9l (x) > 0, (1.1c) 

independent of k. Here / = f(x,t) denotes the density of buyers, g = g(x,t) the density of 
vendors and k the transaction rate. The parameter a is the transaction cost for buyers and 
vendors, a the diffusivity. Note that the densities / and g depend on the parameter k. The 
volume of transactions at a price x is given by: 

V(x,t) = kf{x,t)g(x,t). (1.2) 

Note that we have conservation of the number of buyers and vendors, i.e. 

f(x,t)dx= / fj(x)dx and / g(x,t)dx = / gj(x)dx, 



for all times t > 0. 

The mathematical modelling of (1.1) was inspired by a price formation model presented by 



Lasry and Lions (2007). They considered the same situation described above, but proposed 



a parabolic free boundary problem to model the evolution of the price. The model by Lasry 
& Lions reads: 

r 2 



ft(x,t) 
9t(x,t) 



<j- 



fxx(x, t) + X(t)S(x - p(t) + a) for x < p(t) and f(x, t) = for x > p(t) (1.3a) 
gxx(x, t) + X(t)5(x — pit) — a) for x > p(t) and g(x, t) = for x < p(t), (1.3b) 



where the free boundary p = p(t) denotes the agreed price of the trading good at time t 
and X(t) = —f x (p(t),t) = g x (p(t),t). Note that the difference of buyer and vendor densities 
v = f — g satisfies the following equation 

v t (x, t) = v xx {x, t) + X{t)(5(x - p(t) + a) - S(x - p(t) - a)). (1.4) 

The function X(t) is the transaction rate at time t, which corresponds to the flux of buyers 
and vendors. The Dirac deltas correspond to trading events which take place at the agreed 
price p = p(t), shifted by the transaction cost a. The analysis of system (1.3) was studied 



in a number of papers, see Markowich et al. (|2009|); Chayes et al. (|2009|); Caffarelli et al. 
fl2011a|b| ). 



The main difference between system (1.1) and (1.3) is that the agreed price p = p(t) enters as 
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a free boundary in (1.3). Furthermore the density of buyers and vendors are zero if the price 
is greater or smaller than the agreed price. We shall show that solutions of the Boltzmann 
type equation (1.1) converge to solutions of (1.3) as the transaction rate tends to infinity, i.e. 
k — > oo, thus giving a mathematical justification for (1.3) as a pure formation model. 
This paper is organised as follows: we start with a detailed presentation of the mathematical 
modelling of (1.1) in Section [2j In Section [3] we show that solutions of (1.1) converge to 
solutions of (1.3) as k — > oo. The initial layer problem for (|1.1[ ) is discussed in Section |4j the 
scaling limit k — > oo, a — > and ka = c > in Section [5j This limit corresponds to high 
frequency trading, where computers trade goods on a rapid basis with little or no transaction 
costs involved. The behaviour of both models is illustrated by numerical simulations in Section 

m 



2. Modelling: From agent behaviour to Boltzmann-type equations 

In order to describe price formation at a reasonably simple agent-based level, we consider 
a setup of a large number, say N, of buyers and a large number M of vendors. With the 
exception of discrete events, the price changes are subject to random fluctuations, which we 
model as Brownian motions with diffusivity a. The trading events can be modelled like kinetic 
collisions. After the usual kinetic limit N — > oo, M — > oo we associate the density / = f(x, t) 
to the group of buyers at time t, and the density g = g(x, t) to the group of vendors at time 
t. Here f(x, t) and g(x, t) denote the fraction of the buyers and vendors willing to trade in an 
infinitesimal interval around x at time t. If a buyer with state x meets a vendor with state 
y, they will trade with a certain probability depending on their state. After collision, i.e. 
trade, a buyer becomes a vendor and vice versa. Due to conservation and indistinguishability 
we do not need to take into account the change of roles, but just model a standard collision 
(x,y) — > (x',y') with collision kernel K = K(x,y,x',y'). More precisely, K(x,y,x',y') counts 
the number of trading events per unit time of buyers willing to buy at price x and reselling 
after the trading event at price x' , and the number of vendors willing to sell at price y and 
re-buying after the trading event at price y' . Using the notations for densities as above we 
arrive (with appropriate time scaling) at 

ft(x,t) - ^-f xx (x,t) = / / / K(x',y',x,y)f(x',t)g(y',t) dy dx dy 



2 



a 2 



K(x, y, x', y')f(x, t)g(y, t) dy dx' dy' (2.1a) 




9t(x,t) - —g xx (x,t) = / / / K(x',y',y,x)f(y',t)g(x',t)dxdx'dy' 
^ Jr Jr Jr 

K(y, x, x, y')f(y, t)g(x, t) dx dx' dy' . (2.1b) 

The peculiar aspect in the modelling of the collisions related to the Lasry-Lions approach 
is introduced by the transaction costs. Right after selling at a certain price, it only makes 
sense to re-buy at a price that is at least lower than the previous execution price minus the 
transaction costs, and similar reasoning applies to the buyers. Thus, if we assume symmetric 
transaction costs a and denote by r(x, y) the price at which a buyer with state x and a vendor 
with state y trade, then 

x' = r(x,y)-a, y =r(x,y) + a, (2.2) 
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and hence K is of the form 

K(x, y, x', y) = K (x, y)5(x - r(x, y) + a)S{y' - r(x, y) - a), (2.3) 
where 5 denotes the Dirac 5-distribution. The detailed properties of Kq and r now depend 



on the modelling of the limit order book (cf. Cont and De Larrard (2011); Cont et al. (2010); 



Preis et al. (2006); Bouchaud et al. (2002); Abergel and Jedidi (2011)). Here we make a 



simple assumption that x and y can be interpreted as the bids of buyers and vendors, which 
are tried to be matched exactly. Hence, Kq is centred around x = y, and the simplest choice 
is to assume 



K (x,y) = k5(x - y), 
where the constant k is the trading frequency. We only need to specify 

r(x, x) = x. 

This leads to the collision kernel 

K(x, y, x', y) = 8{x - y)5{x - r(x, y) + a)5{y - r(x, y) 
For smooth test functions (p = (f(x, y, x' , y') we have 

(K,(p) = [ ip(x, x, r(x, x) — a, r(x, x) + a)dx 
ip(x, x, x — a, x + a)dx. 



(2.4) 



Applying this to (2.1) leads to (1.1). For large numbers N and M it appears natural to 
consider the case of large k, since a high potential for transactions is available. One observes 

cfg 

with c 



that the volume of transactions in (1.1) is given by kfg, hence the scaled density p 

yields a density of currently traded prices. If p is well centred, the mean, 



Im f9 dx 

median, or maximum p will give an estimate of the price p(t), which we will also exploit in 
numerical comparison with the Lasry-Lions model. For large k we will see that p concentrates 
to a Dirac 5 distribution centred at p{t) . Note that the mean supply and demand price satisfy 

d 



— I xf(x,t)dx = —a I kf(x,t)g(x,t)dx<0 



xg(x,t)dx = a / kf(x,t)g(x,t)dx>0. 



dt 

This means that the average demand price is decreasing, while the average supply price 
increases as long as fg does not vanish. This shift is clearly inherent in the collision model 
for the trades, since each collision decreases the demand price of an agent and increases the 
supply price of another one. 

We finally mention that the above modelling is reasonably simple to study basic effects, but 
the agent-based interpretation easily allows various generalisations that can make it more re- 
alistic. An obvious first example is replacing the Brownian motion by more general stochastic 
differential equations, which simply changes the differential operator to some other parabolic 
equation, but leaves the collisions unchanged. Even herding effects could be modeled this 
way, e.g. by making the drift of a buyer dependent on the (empirical) density /, e.g. on 
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J" R x/(x, t) dx. The extension to more complicated trading models is inherent in our model 
by specifying Kq and r. A simple alternative model is 

x ~\~ y 

K (x,y) = kl{ y < x y, r(x,y) = (2-5) 

which corresponds to a self-organised trading with using the mean value of the prices when 
y < x. 



3. Limiting equations as k — > oo 



In this section we show that the difference / — g of solutions of the system (1.1) converges 



to the solution of ( |1.4[ ) as the transaction rate k tends to infinity and that, as k — > oo, 
/ = (/ — g) + and g = (/ — g)~ . To emphasise the dependence of / and g on k, we use the 
notation f k = f k (x,t) and g k = g k (x,t) throughout this section. Without loss of generality 

2 

we set = 1 in the rest of the paper. Note that existence and uniqueness of a non- negative 
solutions of (1.1) for k > follows trivially from the estimates stated in this section. For the 
following we assume that // and gi are independent of k and satisfy: 

(A) Let fj, gj > on R and fj,gj G S(R). 

Note that assumption |(^4)| can be weakened, i.e. 

G L X (R) n L°°(R) n iJ x (]R) with sufficiently fast decay as |x| -> oo, 

at the expense of additional technicalities. 

Proposition 3.1. Lei assumption \(A)\ be satisfied. Then system (1.1) /ias unique positive 
solutions f k ,g k G L°°(0, oo; L 1 (M) n L°°(M)) bounded uniformly as k — > oo. Furthermore 

f k g k ^0 inV'(Mx [0,oo]). (3.1) 

Proof. Note that the function u(x,t) = f k (x,t) + g k (x + a, t) satisfies the IVP for the heat 
equation 

ut(x,t) = u xx (x,t) (3.2) 
u(x, 0) = //(x) + + a). 

The heat equation has a unique positive solution u = u(x,t) G L°°(0, oo; L°°(IR)) for initial 
data fi,gi G <S(M), gi{x) > 0, > independent of k. Hence we deduce that the weak 

limits /°°, g°° satisfy 

/°°(-,t) G L+(K) n L^(M) and g°°(-,t) G L+(K) n L~(M) uniformly in t > 0. 

To show the weak convergence of to zero we rewrite the first equation in ( |1.1[ ) as 

i/ t fc (x, i) - i/^(x, t) = t)/(x, i) + / fc (x + a, t)/(x + a, t). 

We define C fc (M) = f k (x,t)g k (x,t), Q k G L°°(0, oo; (R) n L°°(M)). Then in the limit 
k — > oo we deduce that 

— C fc (x, t) + ( k (x + a,t) = o(l) G P'(K x [0, oo)) (^(x, t) = C°°(x + a, t). 

The only periodic function in L 1 (M) is the zero-function, which concludes the proof. □ 
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Next we introduce as in Caffarelli et al. ( 2011a|b ) the functions 

oo oo 

F k (x,t) = ^2f k (x + al,t) and G k (x,t) = ^g k (x-al,t). 



(3.3) 



1=0 



1=0 



Proposition 3.2. Let assumption (A) be satisfied. The functions F and G defined by the 
infinite series (3.3) satisfy F k ,G k 6 L°°(0, oo; L°°(M.)). The series converge locally uniform 
in x and t as k — >■ oo. 

Proof. Since f k (-,t), g k (-,t) £ <S(M) we can estimate the series F k by the following integral 



I(x) 



dy 



o (l + |x + y|) 2 J x (l + \z\) 



-dz for some c > 0. 



l l 



For x > the integral I(x) — (i+x) 
following estimate for the integral 



tends to as x — > oo. In the case x < we obtain the 



I{x) 



1 



:dz 



kl 

< ci + 



-w(i + ^) 2 7o (i + *) 



C2 



- , where c\ , C2 are independent of x. 



+ M 

Note that F k and G k satisfy: 

F k (x,t) = -kf k (x,t)g k (x,t)+F k x (x,t) 
G k (x, t) = -kf k (x, t)g k (x, t) + G k xx (x, t). 



□ 



(3.4a) 
(3.4b) 



Then the difference $ := F k — G k satisfies the heat equation &t = &xx- Note that $ is 
independent of k for fi,gi independent of k. From system (3.4) we deduce that /i k = kf k g k 
converges as k — > oo to a locally bounded positive measure on M x [0, oo). Then the limiting 
functions f°° and g°° satisfy the following system 



^ (x, t) = -v°°(x, t) +^{x + a, t) + f£{x, t) 
5t °° (x, t) = (x, t) + n°° (x -a,t) + (x,t). 



(3.5a) 
(3.5b) 



Note that the difference v = f°° — g°° satisfies a parabolic PDE with a similar structure as 
the price formation model (|1.4|), i.e. 

(3.6) 



v t (x,t) = v xx (x,t) + ii°°(x + a,t)- (i°°(x - a,t). 



3.1. A priori estimates. We derive various a priori estimates for the solutions f k and g k 
of system (1.1), which will we used for the identification of the limiting system as k — > oo. 



Proposition 3.3. The solution of (1.1) satisfies the following a-priori estimate: 

cT 

(f k (x,t) + g k (x,t)) 2 dxdt < const, 



o 



(3.7) 



uniformly as k — ^ oo. 
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Proof. In the following calculations we neglect the dependence of the functions / , g k ,u,F k 
on the variables x and t, and only state their arguments if necessary. We multiply system 
(1.1) with f k and g k respectively and obtain 

\jt J(f k ) 2dx = - k J{f k ?9 k dx + k J f k (x + a,t)f k (x,t)g k (x + a,t) dx- J(f k ) 2 dx, 

\jt J(9 k ) 2 dx = -k J f k (g k f dx + k J f k (x-a,t)g k (x,t)g k (x-a,t) dx - J (g k f dx. 

With changes of variables in the second integral on the right hand side of both equations we 
deduce for the sum of both equations that 

lit! {{fk)2 + {ak)2)dx = ~ k S [fkgk{fk + gk) " fkgk{fk{x - °' t] + gk{x + a> t))]dx 

- f(g k + f k ?dx. 

Next we use that g k (x + a,t) = u(x, t) — f k (x, t) and f k (x — a,t) = u(x — a,t) — g k (x, t) and 
obtain 

1 d 



2dt 



(3.8) 



((/fc) 2 + {gk) 2 )dx = _ k j [fkgk{f k + gk) + {fk) 2 g k + fk {g kf ]dx 

+ J kf k g k (u(x, t) - u(x - a, t))dx - j \g k + f x fdx. 

Multiplication of ( |3.4[ ) by u = u(x,t) and, resp., by u(x + a,t) gives 

J (F k u) t dx - J F k u t dx = J F k u xx dx -k j f k g k udx. 

Note that all integrals are well defined, since u = u(x,t) decays algebraically as |x| —> oo. 
Hence we obtain after integration over the time interval (0, T) that 
rT r 

kf k (x,t)g k (x,t)(u(x,t) + u(x + a,t)) dx dt < const (3.9) 



uniformly as k — > oo. Using (|3.9l) in (|3.8|) concludes the proof. □ 



Multiplying the first equation in system (3.4) by a test function ip G C^°(M x [0, oo)), <p > 0, 
we obtain 

jj kf k g k ipdxdt = jj F k ip t dxdt + / F k ip \J =0 dx + j j F k ip xx dxdt < const. 
If (f(-,t) has compact support on (— 3R, 3R) and tp = 1 in (— 2R, 2R) we deduce that 

rT r 



g (x, t)dxdt < const 

J\x\<R 

uniformly in k. Then the functions f k and g k satisfy 

ft k ,g? £Lj oc (Rx [0,oo)) + L / 2 oc (0,oo;i/- 1 (M)), (3.10) 

and 

f k ,g k GLl c (0,oc;HUR)) (3.11) 

uniformly in k. 
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Strong convergence. Let be an interval on the real line f2 = (—R,R). To show strong 
convergence of (3.1) we use a generalised version of the Aubin-Lions lemma, cf. Lions (2002). 
We define the following spaces 

H 1 ^) CC L 2 {fl) C ff- x (0). 

=B =B =B 1 

For a function v 6 L 1 (0,T; we estimate the corresponding norm by 



(3.12) 



/ mpdx 



o <p<=Hl(n) W\\h1(q,) 



sup 



dx dt < 



T 

sup 

o <p&Hi(n) 

T 



f\u\dx cftlMlz/x^) 



dx dt 



< c 



o 



u(x, t)\dxdt. 



Due to (3.10) and the previous estimate we deduce that f k ,g k e ^(O, T; iT _1 (J2)). Since / fc 
and <7 fc converge weakly to /°° and g°° in L 2 (0, T; L 2 (S7)), we use the general version of the 
Aubin-Lions lemma, cf. Lions (2002), to obtain that 



lim (f k ,g k 

k— >oo 



(r,g°°)^ (L 2 (0,T;L 2 (Sl))f 



(3.13) 



for every bounded interval Hence we conclude strong convergence of the product in 
L^(0,oo;R),i.e. 



k r .k 



lim f K g 

k— >oo 



rg°° = 0. (3.14) 

3.2. The limiting equations as k — > oo. Now we are able to identify the price formation 
model given by Lasry & Lions (1.3) in the limiting case k — > oo. We have already shown 
that the Boltzmann type price formation model has a similar structure in the limit k — > oo, 
see (3.6). Therefore it remains to show that the measures on the right hand side of (3.6) are 
indeed the transaction rates defined by Lasry &; Lions. We shall do this for the initial data 
// and gi, which satisfies the compatibility conditions necessary for (1.3), i.e. 

(5) sup{//(x) > 0} < M{ gi (x) > 0}. 
We find by applying the Lebesgue's dominated convergence theorem that 

oo oo 

F°°(x,t) = J2f°°(x + al,t) and G°°(x,t) = ^g°° (x - al,t). 

1=0 1=0 

Then the difference <3?(x,i) = F°°(x,t) — G°°(x,t) solves the IVP for the heat equation 

$t(x,t) = A$(x,t) (3.15a) 



<f>(x,0) = $ 7 (x) :=Y / f I (x + al)-J2g I (x-al). 



(3.15b) 



1=0 



1=0 



We start with an additional regularity result. 

Proposition 3.4. The solutions f k and g k of (1.1) satisfy f k ,g k £ L 1 (0, oo; Cb(IR)) for 
almost every time t > 0. This implies the continuity of F k and G k for almost every time 
t > 0. 

Proof. The proposition follows from the compact embedding of H 1 ^) into Cfe(IR), therefore 
f k ,g k £ L^O, oo; C b (R)) for almost every time t. □ 
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We now prove the main result of this section. 

Theorem 3.5. Let fj and gj satisfy assumption \(A)\ 

(i) Then the limiting functions f°° and g°° are given by: 

f°°(x,t) = <S> + (x,t) - <£> + (x + a,t), and g°°(x,t) = ®~(x,t) - ®~(x - a,t). 

(ii) Additionally let \(B) hold. Then f°° and g°° satisfy system (1.3). 

The main ingredient of the proof is: 

Lemma 3.6. Let fj and gj satisfy {(A) Then 

(i) F°° = - G°°) + and G°° = (F°° - G°°)~ . 



(ii) Additionally let (B) hold. Then there exists a unique globally defined continuous func- 
tion p = p(t) (the price), which satisfies 

F°°(p(t),t) = G°°(p(t),t) = and 

Proof, (i) Take (x,t) elx (0,oo) such that F°°(-,t), G°°(-,t) are continuous and F°°(x,t) > 
0. Then there exists an I such that f°°(x + al,t) > and therefore g°°(x + al,t) = 0, 
due to ( 3.14| ). The function u(x,t) = f k (x,t) + g k (x + a, t) satisfies the heat equation (3.2) 
and is positive, hence we can deduce that f°°(x — a(l — l),t) > and subsequently that 



g°°(x + a(l — l),t) = 0. By repeating this argument we can show that g c 



al,t) = 



for Z = 0,1,2,..., which implies G°°(x,t) = 0. Analogously we prove G°°(x,t) > implies 
F°°(x,t) = 0. 

(ii) Prom the results in the work by Caffarelli et al. ( 2011a|b| ) we deduce that there exists a 
unique globally defined continuous function p = p(t) (the price) which satisfies 

$(p(t),t) = = F°°(p(t),t) - G°°(p(t),t). (3.16) 

Therefore F°°(p(t),t) = G°°(p(t),t). Assuming that F°°(p(t),t) > leads to a contradiction 
by proceeding similarly to the proof of (i). □ 



This allows us to identify the limiting functions f°° and g°° in the case of (ii) of Theorem 3.5 
and conclude, as in |Caffarelli etH] fl2011a|b|) that f°° an d g°° sati sfy system JL3| ). 
We remark that the theory developed in Caffarelli et al. ( 2011a|b ) applies strictly speaking 
only if sup{//(x) > 0} = mi{gi(x) > 0}, however an extension for sup{//(x) > 0} < 
inf{<7/(x) > 0} is straight forward. 

Remark 3.7. In the case (ii) the transaction volume in the limit k — > oo is 

fi°°(x,t) = X(t)5{x - p(t)) (3.17) 
and the limited density of traded prices is p°°(x,t) = 5(x—p(t)). To see this we compare (|3.6|) 



to (1.4) and obtain: 

.oo 



p,°°(x + a,t)- fi°°(x -a,t) = X(t)(5(x - p(t) + a) - S(x - p(t) - a). 

This implies that fjP°(x,t) = X(t)5(x—p(t))+A(t,x), where A is a nonnegative and 2a-periodic 
function. Now we expand A into its Fourier series and consider a single harmonic term in 
the series, given by ai{t)e mx a ,/ £ Z. When taking this term as inhomogeneity in the heat 
equation 

z t = z xx - ai(t)e inx a 
z(x,t = 0) = 0, 
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we easily compute the solution: 

z(x,t) = -e mxl a f e-^ 2 ^ ai {s)ds. 
Jo 



Passing to the limit k to infinity in (3. 4) (a) gives the heat equation for F°° with — as 



inhomogeneity. From the proof of Proposition \3.%\ it becomes clear that F°° tends to for 



x — > oo (for every fixed t) and thus F°° does not admit x-periodic modes. We conclude (3.17). 



If assumption (B)\ does not hold, then it can be shown that the local transaction rate 



satisfies in the limit k — > oo: 

H°°{x,t)= \^{pj{t),t)\8(x- Pj (t)), 
je.J(t) 

where J(t) denotes the index set (finite or countably finite) such that {pj{-,t) \ j G J(t)} is 
the set of zeros of $(•,£). 

4. The initial layer problem - preparation of the initial data 
In this section we discuss the behaviour of the initial layer on the fast time scale. This 



initial layer occurs if (f°°(x, t = 0),g°°(x, t = 0)) as computed in Theorem 3.5 (i) differs from 
(fi,9i)- 

Let e = then system reads 

sf*(x, t) = -f £ (x, t)g £ (x, t) + f £ (x + a, t)g e {x + a,t) + ef £ xx (x, t) 
£9 e t(x, t) = ~f{x, t)g e (x, t) + f{x - a, t)g E (x - a,t) + eg £ xx (x, t). 

Let t = - denote the fast time scale and the corresponding fast-scale dependent variables by 

a £ (x, t) := f £ (x, t) and f3 £ (x, r) := g £ (x, t). 

Then a £ = a £ (x,r) and (3 £ = f3 £ (x,r) satisfy 

a £ (x, t) = —a £ (x, t)(3 £ (x, t) + (3 £ {x + a, r)a £ (x + a, r) + ea xx (x, r) (4.1a) 
/3 £ (x, t) = -(3 £ {x, r)a £ (x, r) + f3 £ (x - a, r)a £ (x - a, r) + £^ x (x, t). (4.1b) 

with initial data 

a £ (x, 0) = fi(x) and /3 £ (x, 0) = g r (x). (4.1c) 

Proposition 4.1. In the limit e — > the fast scale variables a £ = a £ (x, r) and (3 £ = f3 £ (x, r) 
converge to 

a e ^ «0 ; p e ^ p0 in L 2 oc(0)OO;L 2 oc(M)) . 

The limits a = a°(x,r) and /3° = f3°(x,r) satisfy the ODE system 

oP t (x, t) = -a°(x, t)/3°(x, t) + a°(x + a, r)/3°(x + a, r) (4.2a) 
^(x, r ) = -a°(x, t)P°(x, t) + a°(x - a, r)/3°(x - a, t). (4.2b) 



Furthermore the system (4.2) with IC (4.1c) has a unique continuous space-time solution. 
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Proof. Note that v £ (x,t) = a £ (x,r) + f3 £ (x + a, r) satisfies the heat equation 

v £ t (x,t) = ev £ xx (x,r) 
v £ (x,0) = fi{x) + g I {x + a). 
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(4.3a) 
(4.3b) 



Solutions of (4.3a) with initial datum fi(x), gi(x) G S(R) decay algebraically fast and are 
smooth, and we deduce t hat a £ , f3 £ G L°°(0, oo; L X (R) n L°°(R)) uniformly as e -> 0. Next 
we differentiate equation (4.1a) with respect to a;, multiply with a| and integrate over R to 
obtain 

1 d 



(o4) 2 cte < K / (a^cte + if / {figfdx - e I (a £ xx )' z dx. 



2dr 

The same holds for ( 4.1b| ), i.e 

2d7 



,.e\2 



dx<K [ (a £ x ) 2 dx + K [ (f3 £ x ) 2 dx 
Jr Jr 

Integration over (0, T) gives, for arbitrary T > 0: 

[(a £ x ) 2 (x,T) + ^ £ x ) 2 (x,r)]dx<K 1 (T) 



e / (f3 £ xx ) 2 dx. 



for r G (0, T\. Similar calculations, i.e. differentiating equations (4.1a) and (4.1b) with respect 
to r and multiplication by a% and (3 £ respectively, lead to 



[(a £ T r(x,r) + (f3 £ r(x,r)]<K 2 (T), 

for r G (0, T). Therefore a £ and /3 £ are bounded in Hj oc (R x (0,oo)) uniformly as e — )• 0. 
Due to the compact embedding of H} oc into Lf oc we can deduce that, after extraction of a 
subsequence 

of a , (3 £ -»• /3° in L 2 OC (0, oo; L 2 oc (R)). 

The existence of a unique continuous space-time solution of (4.2), given by (a°,j3°), is imme- 
diate. □ 

Next we study the behaviour of ct° = a°(x, r) and j3° = j3 (x, r) as r — > oo, which corresponds 
to the "end of the initial layer" . 



Proposition 4.2. (a) Solutions a = a (x,t) and /3° = j3°(x,r) of (4.2) converge in the 
limit t — > oo 

lim q°(x,t) = a°°(x), lim /3°(x, r) = /3°°(x) and a 00 (x)/3 00 (x) = onR. 

T—^OO T—>00 

(b) If f I (x)g I (x) = onR then a°(x,t) = fi(x) and (3°(x,t) = g^x) on R x [0,oo). 

Proof Set w° = a°(x, t)(3°(x, t) and compute 

w°(x, t) = a°(x, t)w°(x - a, t) - [a°(x, r) + (3°{x, t)]w°{x, t) + /3°(x, t)w°(x + a, r) 

w°(x, 0) = a°(x, 0)/3°(x, 0) = fi{x) gi {x). 

If fi(x)gi(x) = on R, then u>°(x,r) = a°(x, r)/3°(x, t) = 0. Hence (4.2) gives a°(x,r) = 
fi(x) and /3°(x, r) = 57 (x) on R x [0, oo). 
We observe that 

a°(x, r) + f3°(x + o, r) = //(a;) + + a) 
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and 

oo oo oo oo 

a°(x + al, t) — ^2 P°( x ~ T ) = 5^ fl( x + a U T ) — SlO^ ~~ a h T )- 

1=0 1=0 1=0 1=0 

v v ' s v ' 

=:A°(x,t) =:B°(x,t) 

Note that a (•, r), /3 (•, r) G thus both sequences A (-,r) and B°(-,t) converge uni- 

formly locally in x. Moreover both sequences are decreasing in time, since 

A° T = -w° < and B° T = -w° < 0. (4.4) 

Therefore for every i£R fixed, we have 

A°(x,t) 4 A°°(x) and B°(x,t) | B°°(x) as r -> oo. 



Integration of (4.4) in time gives 



A°(x,r) - A°(x,0) = / A° T (x,s)ds = - w°(x,s)ds. 

Jo Jo 

Therefore we conclude from the existence of the limits A°°(x) and B°°(x) for r — )• oo that 

w°(-,t) £ L^_(0, oo) for every 
Then a^(x, r) = — w°(x, r) +ui°(i + a, t) and integration over the interval [0, r] gives 

rr rr 

a°(x, t) — fj(x) = — / u;°(x,s)ds + / w°(x + a, s)ds. 



We know that u>°(x, •) G £+(0, oo) and therefore we conclude for r — > oo that 

lim a°(x,r) = a°°(x), lim /3°(x,r) = /3°°(x) and a 00 (x)/3 00 (x) = w°°(x). 



,,0/ V T^OO. 



Thus a^(x,r) >• — w°°(x) + w°°(x + a) = 0, and we conclude that a°° (x) (3°° (x) = on 

E. ' □ 

The properties of a°° and can be can be summarised as follows: 
(PI) a°°(x) + /3°°(x + a) = //(x) + ^/(x + a). 

oo oo oo oo 

(P2) a °°( x + al ) - J2 /9°°(z - a/) = £ //(x + a/) - £ 5/ (a; - a/). 

i=0 «=0 1=0 1=0 

V v ' V v ' 

=:A°°(x) =:B°°(x) 

(P3) a 00 (x)/3 00 (x) = 0. 
(P4) a°°(x) > 0, p°°(x) > on E. 
Furthermore we set C°°(x) = A°°(x) - B°°(x). 



Theorem 4.3. Let fi,gi satisfy (A) and h(x) := fi(x) + gj(x + a), x G K satisfy 

(i) if for some x\ G M : h(x\) > and i/ i/iere is an li G N wrei/i hj(xi — all) = 0, then 

/i/(xi - a(/i + 1)) = 0, 
(wj if for some x 2 G R : h(x 2 ) > and i/ i/iere is en I2 6 N wii/i hi{x 2 + a/2) = 0, then 
h I (x 2 + a(l 2 + l)) = 0. 

Sei $/(x) := X]£o/^( x + a ~~ YH^=o9i( x ~ a 0- I7ien i/ie limiting functions a°°(x) = 
limT-^oo a°(x, r) and /3°°(x) = lim^oo /3°(x, t) can 6e identified as: 

a°°(x) = $+(x) - $|(x + a) and /3°°(x) = $7(3) - $7(x - a). (4.5) 
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Proof. The proof proceeds in analogy to the proof of Proposition 3.6 and is split into two 
parts. First we show that the positivity of either a 00 or (3°° at x G M implies that the other 
variable is zero at x and x =L al as well. In the second step we identify the limiting functions. 

(a) Let x G R be such that a°°(xo) > 0. Then /3°°(xo) = for all I G N, because of 
property (P3). Next we show by induction that /3°°(xo — al) = 0, V7 G N. Therefore 
let for h G N hold /3°°(x - al) = for all I G {0, 1, 2, . . . , h - 1}. Then we have 

a°°(x - ah) + /3°°(x - a(h - 1)) = /i/(x - ah) (4.6) 

V v ' 

=0 

and 

a°°(x ) + /3°°(x + a) = h^xo) > 0. 
We distinguish between two cases: 



Case 1: If hi(xo — ah) > 0, then a°°(xo — ah) > (because of (4.6)). Hence we deduce, 

using property (P3), that /3°°(xo — ah) = 0. 
Case 2: If Iii(xq — ah) = 0, then by assumption (i) of the theorem we have h{xo — a{h + 
1)) = 0. Property (PI) gives a°°(x - a(h + 1)) + P°°{xq - ah) = hi(x - a(h + 
1)) = and we conclude that (3°°(xq — ah) = 0. 
(b) Let xq G M with (3°°(xq) > 0. From property (P3) we conclude that a°°(xo) = 0. 
Again we use an induction argument to show that a°°(xo + al) = for all / G N. 
Therefore for I2 G N let a°°(xo + al) = for all / G {0, 1, . . . , h — 1}- In this case 

a°°(x + a{l 2 - 1)) +/3°°(x + al 2 ) = h^xo + a(l 2 - 1)) (4.7) 



and 



>(x -a)+ f3°°(x ) = h T (x - a) > 0. 



Again we consider the two different cases. 



Case 1: If hi{x + a{l 2 - 1)) > 0, then /3°°(x + al 2 ) > (because of ( |4~7| ) ) and therefore 

a°°(xQ + al 2 ) = follows. 
Case 2: If hj(x + a(l 2 — 1)) = 0, then by assumption (ii) of the theorem we know that 

hj{xQ + al 2 ) = and therefore a 00 (x + al 2 )+f3 00 (xo + a(l 2 + l)) = hj(x + al 2 ) = 0. 

Thus a°°(x Q + al 2 ) = 0. 

Now we identify the limiting functions a°° = a°°(x) and f3°° = (3°°{x). We choose an xq G M. 
such that ^°°(xo) > 0. Then there exists l 2 G N U {0}, such that a°°(a;o + al 2 ) > 0. From 
(a) we deduce /3°°((x + al 2 ) - al) = /3°°(x - a(l - l 2 )) = for all I G N. Therefore 
^(xo) = 0. Analogously B°°{x 1 ) > implies A°°(x{) = 0. Therefore A°°(x) = C°°(x)+ 
and B°°(x) = C°°(x)- , i.e. 

A°°(x) = $|(x) and B°°(x) = §J(x). 

We conclude that 

a°°(x) = $+(x) - $f(x + a) and p°°(x) = $J(z) - $j(x - dp 



Under the assumption of Theorem |4.3| on fj and gj we computed the unique solution 
candidate satisfying properties (P1)-(P4). Since existence of a solution of (P1)-(P4) follows 
already from the r — > 00 argument, we conclude that we have indeed calculated a 00 ,/? 00 . 
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Figure 1. The initial datum depicted in figures (a)-(c) satisfies the assump- 



tions of Theorem 4.3 This is not the case for the initial datum depicted in 
Figure (d). 



Note that Theorem 4.3 does not cover initial functions // and gi £ C^°(IR), whose supports 
are separated by a distance larger than a (see Figure [T] (d)). If the distance between the 



supports of // and gj is smaller than a or if they intersect, Theorem 4.3 is applicable. 



Remark 4.4. Note that v 6 (x,t) = a £ (x,r) + (3 s (x + a, r) satisfies the heat equation (4.3), 
thus 

lim (a e ,/3 e ) = in L°°(R), 

T — ► OO 

and, actually: 

< a £ (x,r) + P £ (x,t) < HMInnfl Vx G M > Vr > °- 

Therefore the limits r — > oo and e — > do not commute. Note that the initial layer problem 
(4.1 ) does in all generality not provide well-prepared initial data satisfying assumption \(B)\ We 



conjecture that small diffusion helps in his process, i.e. that a £ , j3 £ are close to well-prepared 
on large time scales r = r(e) with er(e) small (see Example 3 in Section^. 

5. Motivation by fast trading with small transaction rates 
In this section we consider the scaling limit k — > oo, a — > with ka = c = const. We rewrite 



system ( 1.1 ) as 



/«(*, t) = c to)(* + °.*) -(*?)(*■') + f xx{x , t) (5.1 a) 

a 

i t\ (fg)( x ~a,t)- (fg)(x,t) 

g t {x,t) = c \-9xx{x,t). (5.1b) 

a 

Formally it is clear that the limiting system as k — > oo and a — > is 

t) = c(f°g°Ux, t) + f^(x, t) (5.2a) 

g° t {x,t) = -c(f°g Ux,t)+g° xx (x,t). (5.2b) 
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Note that u (x, t) = f (x,t) + g°(x,t) satisfies again the heat equation u^(x,t) = u xx (x,t) 
with initial datum u°(x, 0) = fi(x) + gi(x). The rigorous statement for this limit is stated in 
the following theorem: 

Theorem 5.1. Let k > ; a > ; ka = c and fi,gi € Then the weak limits f°, g° 

satisfy system (5.2) subject to /°(x,0) = fi(x) and g°(x,0) = gi(x). 



Proof. We reiterate that u(x, t) = f(x, t)+g(x+a, t) satisfies the heat equation (3.2), therefore 
we can deduce that 

IMIl°°(o,oo;]R) < K uniformly as k — > oo, a — > with ka = c. 

Multiplication of system (5.1) with / and g respectively gives 

f(x + a, t)g(x + a, t)f(x, t) - f 2 (x, t)g(x, t) 



1 d 
2dt 
1 d 
2dt 



f dx = c 



g dx = c 



dx 



f(x + a, t)g(x + a, t)g(x, t) - f(x, t)g 2 (x, t) 



dx — i g x dx. 



Now we look at the first term on the right hand side of the equation for / and deduce 

f(x + a, t)g(x + a, t)f(x, t) - f 2 (x, t)g(x, t) 
ax 

: a 

ff f . ,w/ g{x + a,t) - g(x,t) f(x + a,t)-f(x,t) 
[f(x + a, t)f(x, t) h f(x, t)g{x, t) }dx 



< 



L°° 



l/(M)|| 



g(x + a,t) -g(x,a) 



dx + \\g\\L°° / \f{x,t)\\ 



J(x + a,t)-f(x) 



\dx 



<i||/IU-i 
+ fll/IU- 



1 If 

f 2 (x,t)dx + -\\g\\ L ^- / f 2 (x,t)dx 

g(x + a,t) -g(x,t) 2 fi 

dx+-\\g\\L° 

a I 



f(x + a,t)-f(x,t). 



'dx. 



The difference quotient for / and g can be estimated using the Fourier transform, i.e. 



f(x + a,t) - f(x,t) 2 



dx= / 



l di. 



is bounded by |£| we obtain that 
f(x + a,t) - f(x,t) 2 



dx < L \f x (x, t)\ dx. 



Similar computations hold for the second equation in g. Now we choose ji 
then there exists a a > such that 



2(ll/lli 



2dt 



(f(x,t) 2 + g(x,t) 2 )dx + a / (/s(M) 2 + 2*M) 2 )<te< M. 



Thus 



f,g e Lf oc (0, oo; iZ^QR)) uniformly for a ->• 0. 



By the Aubin Lions lemma, Showalter ( 1997), we conclude that / — > f° and g — > g° as a — > 0. 
Thus we can pass to the limit in system (5.1), making the formal limit rigorous. □ 
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Remark 5.2. Consider the price formation FB problem (1.4) subject to the initial condition 
v(x, 0) = fi(x) — gi{x), where fj and gj satisfy \(Aj\ and (B), Being interested in the limit as 
the transaction fee a tends to zero we denote v = v a , f = f a and g = g a ■ The initial data fi 
and gi are assumed to be independent of a. It is an easy exercise to show that the weak limits 



> v°, f a ^> f° and g a 



a->0 



d 



F\f° = (v°) 



> g° exist and that 
d 



:(F°) 



dx dx 
holds, where F° satisfies the IVP for the heat equation 



,9 



,.0\- 



dx 



0\- 



F° 



F°(x,0) 



x £ 

oo 



,t > 



fl(x + y)dy 



o 



9i{x - y)dy. 



Using the methods of Caffarelli et al. (2011afb ) we conclude the existence of a unique contin- 
uous function p° = p u (i) (the price of the limiting system) such that 



F°(p°(t),t) 



0. 



The parabolic Hopf lemma implies ^F°(p°(t),t) < 0. Therefore the function v°(-,t) has a 
jump discontinuity at x = p°(t), with 

v°(x,t) < 0. 



lim v°(x,t) 
x-*p°(t)+ 



lim 

V(t)- 



We conclude that v a (x,t) has an internal layer of rapid transition at p = p a (t) for a small. 
Moreover it follows that the limits k — > oo (from the Boltzmann type model to the Lasry Lions 
FBP) and a — > (from the Lasry Lions FBP to the limit as discussed in the remark) consec- 
utively do not give the same result as k — > oo ; a — > and ka = c > in the Boltzmann type 
model. 

We conclude that for markets with many transactions with small fees, the Boltzmann price 
formation model is very sensitive to the relative sizes of the transaction frequencies and trans- 
action fees. We refer to the next section for numerical illustrations. 



6. Numerical simulations 



In this section we illustrate the behaviour of the Boltzmann type price formation model (1.1) 



on a bounded domain ft. We supplement system ( 1.1 ) with homogeneous Neumann boundary 



conditions, hence the total mass of buyers and vendors is conserved over time. We simulate 



the behaviour of solutions of (1.1) for different values of k and compare them to the solu- 



tions of (1.3). Furthermore we perform numerical simulations in the case of not well prepared 



initial data and study the evolution of the initial layer on the fast time scale, see also Section|4j 

All simulations are performed on the interval Q, = (0, 1), where 1 corresponds to the scaled 
maximum price of the traded good. The domain Q is split into iV intervals of size h. We 
denote the discrete grid points by xi and the time steps by At. Then the discrete solutions 
/" and g™ correspond to the functions / and g at time t n = nAt and Xj = jh. 



The simulations of the free boundary problem (1.3) are based its formulation as an IVP for 
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the heat equation, i.e. 



Vt(x,t) = V xx (x,t) for all x 6 

V x (0,t) = V x (a,t) and V x (l,t) = V x (l-a,t), 



(6.1a) 
(6.1b) 



with v = v(x, t) being a solution to ( |1.4[ ) and V given by 

V(x,t) = 



Y,n=o v + {x + al,t) 



Caffarelli et al. (2011a). Note that ( 6.1b|) c orresponds to the transformed homoge- 



as in 

neous Neumann boundary conditions. System (6.1) is solved using an implicit in time finite 
difference method. The solutions / and g to ( |1.3[ ) can be calculated from 

v(x,t) = V(x,t) - V + {x + a,t) + V~(x - a,t). 

The price p = p(t) corresponds to the zero level set of V = V(x, t). 

System (1.1) is solved using a semi-implicit in time discretization, i.e. the diffusion is dis- 
cretised implicitly in time, the "collision" terms explicitly. This results in the following 
discretization: 



At 



\ 



At 



(6.2a) 
(6.2b) 



Thus the solution of (6.2) corresponds to solving the linear system 

7 nN 



(I + At A) 



9' 



(/ + AtB) 



n-l 



n-1 



where the matrices I, A, and B are to the identity, the discrete Laplacian and the discrete 
"collision" term respectively. 



Example 1: Comparison of the models for large trading rates k. We compare the 

fj(x) and gi = 
which satisfies the compatibility condition 



behaviour of solutions of (1.1) with (1.3) using the same initial data // 
gi(x). Therefore we choose an initial datum 



W = -§£(P(0)>0) = g (p(0),0), namely 



-10x + 6 



if x < 0.5 

if 0.5 < x < 0.6 

otherwise 



9i{x) 



10x-6 ifx>0.6 
otherwise. 



Since (1.1) converges to (1.3) as k — > oo, we set k = 10 . The other parameters are set to: 

and h = 0.002. Figure 2] shows the evolution of the price and 



-6 



a = 10h, a = 1, At = \ = W- 

the distribution of buyers and vendors at time t = 1 for both models. We observe that the 
buyer-vendor distribution as well as the evolution of the price agree very well. 



IS 
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(a) Buyer and vendor distribution at time 
t = 1 



(b) Evolution of the price p(t) in time 



Figure 2. Comparison of the price formation models for conforming initial data 



Example 2: Behaviour of solutions of (1.1) for different trading rates. Next we 



illustrate how solutions of (1.1) change for different values of k. The initial datum is set to 



//(*) 



9i (x) 



15(x-0.3)(0.5-x) 




15(0.55 - x)(x 




0. 



if 0.3 < x < 0.5 
otherwise 

if 0.55 < x < 0.8 
otherwise. 



The simulations were performed with the following parameters: h = 10 3 ,r 



k 



10 5 and 10 6 , a = lO/i. The final distribution of buyers and vendors as well as the evolution 
of the price is depicted in Figure |3j We observe that the value c = /(x, 0.5) = /(x,0.5) 
decreases if k increases, also the evolution of the price towards its stationary value is slower. 



5 as in (1.3) 



The support of the trading distribution (1.2) decreases for larger k and converges to a Dirac 



Example 3: Not well prepared initial datum on the fast time scale. In this exam- 
ple we study the behaviour of solutions in the case of a not well prepared (in the sense of 



assumption (B) ) initial data fj and gj on the fast time scale. Hence we simulate system (4.1 ) 



with an initial datum, where the density of vendors and buyers are switched, i.e. the initial 
data are not well prepared in the sense of assumption \(B) We set 

(l5(x - 0.65)(0.95 - x) if 0.65 < x < 0.95 
I otherwise 



9l(x) 



12(x - 0.25)(0.5 - x) if0.25<x<0.5 







otherwise. 



We choose a spatial resolution of h = 2 x 10 3 , the time steps are At = 2 x 10 



the 



transaction rate k = 5 x 10 and the transaction cost a = lOh. The evolution of sup f(x) > h 
and the inf g(x) > h as well as the final distribution of buyers and vendors is depicted in Figure 
[4} We observe that even though the initial data is not well prepared a £ and (3 £ converge to a 
well prepared solutions in the sense of assumption (B)\ 
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(a) Buyer and vendor distribution at time t = (b) Evolution of sup f(x) and inf g(x) in time 

(dashed line) and t = 1 (solid line) 



Figure 4. Non- conforming initial data on the fast time scale 



Example 4: Scaling limit. Finally we consider the scaling limit discussed in Section [5] We 
choose a large domain Q = [0, 20] and an initial datum of the form 



//(*) 



20 



if 10 < x > 11 
otherwise, 



1 if 9 < x < 9.5 

-2x + 20 if 9.5 < x < 10 
otherwise 

to avoid any interference with boundary conditions. We observe the evolution of the price for 
the following set of parameters for (1.1) 

a = h 



2 x 10' 5 ,k 



5 x 10 4 ,ct = 1 and r = 2 x 10~ 



This parameter set corresponds to the case where ak = 0(1), i.e. the scaling limit in Section 
[5} Figure [5] shows the solutions of (1.1) and (1.3) at time t = 1 and the evolution of the 
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Figure 5. Non- conforming initial data on the fast time scale 

price p = pit)- We observe that the consecutive limits k — > oo and a — > in Lasry&Lions 
do not give the same result as the scaling limit discussed in Section [5} Also the layer of fast 
transition in the Lasry-Lions solutions / and g is clearly visible in Figure 4a, as mentioned 
in Remark 2 of Section [5j 

Acknowledgements 

LC was partially supported by a grant of the DMS division of the NSF. PAM expresses his 
gratitude to the Humboldt Foundation for awarding the Humboldt Research Award to him, 
which allowed him to spend time with Martin Burger's research group in Munster, where 
this research was initiated. PAM also acknowledges support from the Paris Foundation of 
Mathematics. MTW acknowledges support from the Austrian Science Foundation FWF via 
the Hertha-Firnberg project T456-N23. We thank Bertram During (University of Sussex) for 
the useful hints to literature. 

References 

Abergel, F. and Jedidi, A. (2011) A mathematical approach to order book modeling. Social 

Science Research Network Working Paper Series . 
Bouchaud, J. P., Mezard, M. and Potters, M. (2002) Statistical properties of stock order books: 

empirical results and models. Quantitative Finance 2, 251. 
Caffarelli, L.A., Markowich, P. A. and Pietschmann, J.F. (2011a) On a price formation free 

boundary model by Lasry and Lions. C. R. Math. Acad. Sci. Paris 349(11-12), 621-624. 
Caffarelli, L.A., Markowich, P.A. and Wolfram, M.T. (2011b) On a price formation free 

boundary model by Lasry and Lions: the Neumann problem. C. R. Math. Acad. Sci. Paris 

349(15-16), 841-844. 

Chayes, L., Gonzalez, M.d.M., Gualdani, M.P. and Kim, I. (2009) Global existence and 
uniqueness of solutions to a model of price formation. SIAM J. Math. Anal. 41(5), 2107- 
2135. 

Cont, R. and De Larrard, A. (2011) Price dynamics in a Markovian limit order market. ArXiv 
eprint arxiv:1104-4596 . 



ON A BOLTZMANN TYPE PRICE FORMATION MODEL 



21 



Cont, R., Stoikov, S. and Talreja, R. (2010) A stochastic model for order book dynamics. 
Oper. Res. 58(3), 549-563. 

Lasry, J.M. and Lions, P.L. (2007) Mean field games. Jpn. J. Math. 2(1), 229-260. 

Lions, J. (2002) Quelques methodes de resolution des problemes aux limites non lineaires. Les 
Cours de reference. Dunod. 

Markowich, P.A., Matevosyan, N., Pietschmann, J.F. and Wolfram, M.T. (2009) On a para- 
bolic free boundary equation modeling price formation. Math. Models Methods Appl. Sci. 
19(10), 1929-1957. 

O'Hara, M. (1998) Market Micro structure Theory. Blackwell Business. John Wiley & Sons. 

Preis, T., Golke, S., Paul, W. and Schneider, J.J. (2006) Multi-agent-based order book model 
of financial markets. EPL (Europhysics Letters) 75(3), 510. 

Showalter, R.E. (1997) Monotone operators in Banach space and nonlinear partial differen- 
tial equations, Mathematical Surveys and Monographs, vol. 49. Providence, RI: American 
Mathematical Society. 



